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The aim of this study was to investigate the possible association in the interleukin-6 (IL-6) gene with 
Rheumatoid arthritis (RA) in Chinese Han population from Shandong Province. Target regions of IL-6 gene 
were amplified by polymerase chain reaction (PGR) and genotyped. A logistic regression analysis was 
performed to detect potential associations in our case-control sample, the odd ratio(OR) and 95% 
confidence intervals(CIs) were calculated. Furthermore, we systematically tracked all the published studies 
in the field and performed a meta-analysis for the single nucleotide polymorphisms (SNPs) under study. 256 
RA patients and 331 healthy controls were recruited into the case-control study. We found allele frequencies 
of rsl800795, rsl800797 and rsl474347 in RA patients differ from control subjects (P - 0.016, 0.024, 0.020, 
respectively). Significant difference was observed in haplotype frequencies of GCCGCT between RA 
patients and controls (P - 0.0001, OR - 4.066, 95%CI - 1.891 ~ 8.746), while GGCGCT frequencies was 
found lower in RA than controls (P=0.006, OR - 0.669, 95%CI - 0.501 ~ 0.894). The resuhs of the 
meta-analysis showed association polymorphism within the IL-6 promoter with RA. These findings suggest 
that rare IL-6 gene polymorphisms may associate with RA susceptibUity in Han Chinese populations; 
however further studies are needed to assess the vaUdity of the association of IL-6 with RA. 

Rheumatoid arthritis (RA) is a common chronic autoimmune inflammatory joint disease, characterized by 
erosive symmetrical polyarthritis and systemic extra-articular manifestations, affecting approximately 
0.5% to 1% of the adult population worldwide'. Although the etiology and pathogenesis of RA is not clearly 
defined, a combination of risk alleles from numerous genes contributed to RA susceptibility has been found^. 
Among all these genes, HLA-DR loci was most thoroughly studied yet HLA-DR accounts only for approximately 
one-third of the genetic susceptibility to RA. Other loci also implicated in the predisposition to RA, such as 
CTLA4,PADI4,PTPN22,TNF-cx,IL-ip and STAT4 et al account for relatively small additional contribution RA 
genetic susceptibility^. 

IL-6, one of the most well-studied pro-inflammatory cytokines, is a major mediator of the acute phase response 
and was increased in the synovial fluid and tissue of RA patients. Blockade of IL-6 activity with a soluble anti-IL-6 
molecule tocilizumab has been shown to reduce disease activity and radiological progression in I^^. Data has 
shown IL-6 expression might be in part genetically modulated by polymorphisms located at positions rsl800795, 
rsl800796 and rsl800797 (also known as IL-6 -174G/C, -572G/C, -597G/A) in the promoter region of IL-6''-l IL- 
6 promoter polymorphisms have been associated with susceptibility to RA, but with conflicting results in different 
populations'". 

In this study we sought to assess potential susceptibility of variants in promoter and intron area of IL-6 with 
RA. We designed a case-control study to evaluate these associations in a Han Chinese population. Furthermore, 
we performed a systematic review based on published data on various ancestral groups and performed a meta- 
analysis including our results. 
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Results 

The six SNPs of IL-6. Briefly, two hundred and fifty-six Han Chinese 
patients with RA [F/M: 196/60, mean age 50.26 ± 12.86 years, 
median disease duration 3(0.8 ~ 8.0) years] and 331 healthy 
controls (F/M: 239/92, mean age 48.08 ± 13.92 years) from the 
same geographic and ethnic background were included in this 
study(Supplementary Table 1). No difference in age and sex 
distributions were detectable between the two groups (P > 0.05). 

Genotype and allelic frequencies of the SNPs rsl800795, 
rsl800796, rsl800797, rsl524107, rs2069840 and rsl474347 of IL- 
6 are shown in Table 1. There was a statistically significant asso- 
ciation observed for rsl800795, rsl800797 and rsl474347. 
Specifically the ORs(95% CIs) were OR = 4.823 (95%CI = 1.338- 
17.379, P = 0.016), OR = 5.905 (95% CI = 1.270-27.448, P = 0.024) 
and OR = 3.423 (95% CI = 1.212-9.665, P = 0.020) for rsl800795, 
rsl800797 and rsl474347, respectively. Interestingly the minor allele 
frequencies (MAF) of these variants in the control population of our 
study were different from the reported European groups. In addition, 
no relationship was observed when we analyzed IL-6 genotypes dis- 
tributions according to age, gender, mean disease duration, ESR, 
CRP, DAS28 and serum level of RF, aCCP, or aGPI 
(Supplementary Table 2, Supplementary Table 3). 

Meta-analysis of IL-6 SNPs. There are four major haplotype in IL-6 
polymorphism: GCTGCT, GGCGCT, GGCGGT, GCCGCT. 
Significant difference was observed in haplotype frequencies of 
GCCGCT between RA patients and controls (P = 0.0001, OR = 
4.066, 95%CI = 1.891 ~ 8.746). Haplotype GGCGCT frequencies 
was found lower in RA than controls (P = 0.006, OR = 0.669, 95%CI 
= 0.501 ~ 0.894). (Table 3). 

Our systematic search identified greater than 4 eligible articles 
studying the potential association between a variant and rheumatoid 



arthritis for rsl800795 and rsl800796 only. Fourteen studies (13 for 
rsl800795 and 6 for 1800796) including our case-control study accu- 
mulating a total of 2791 RA patients and 3100 control subjects were 
available in the literature. Information on the primary data is shown 
in Table 2. Specifically the MAFs for rsl800795, rsl800797 and 
rsl474347 in our population were 0.5%, 0.3%, 0.8% compared to 
46.5%, 47.3% and 46.0% in Caucasians. Given the difference in the 
MAF between Han Chinese and Europeans a sensitivity analysis by 
ethnic group was performed. 

In a sample-size of 584/703 case/controls, Asians meta-analysis 
revealed an association between RA and the rsl800795 C allele {OR 
= 9.75, 95% CI = 4.99-19.06, P < 0.00001) with no observed het- 
erogeneity {P = 0) (Figure 1). The p-value was 0.8028 from egger's 
test with asymmetry in funnel plot as alternative hypothesis. 

However no difference was observed in Europeans {OR = 1.00, 
95% CI = 0.84-1.19, P = 1.00) with high heterogeneity {P = 68%) 
(Supplementary Figure 1). Furthermore, no association with RA was 
found for rsl800796 C allele {OR = 1.08, 95% CI = 0.54-2.16, P = 
0.82) in Asian populations with notable high heterogeneity {P = 
92%) (Supplementary Figure 2) and no studies were available in 
European populations. 

Discussion 

Our study and the conducted meta-analysis supports evidence that a 
rare polymorphism of IL-6 may play a major role in the susceptibility 
to RA, at least in the Chinese population. The association was nom- 
inally statistically significant in our study with a strong association 
found when all studies were synthesized and with minimal observed 
heterogeneity. 

The summary effect for this variant is null in European popula- 
tions but we have shown that a different LD pattern may be observed 
between these two ethnic groups. A previous meta-analysis of IL-6 
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that attempted to summarize the effect sizes of all ethnicities^ should 
be interpreted with caution given this diverse LD pattern. Indeed, it 
has been shown that differential ancestral effects can be anticipated 
and genomic risk markers may need separate further evaluation in 
different ancestry groups. 

Polymorphism rs 1800795 in IL-6 consists of a G to C substitution, 
which has been reported to be associated with systemic onset of 
juvenile idiopathic arthritis (SoJIA) and to affect levels of IL-6''. 
Serum levels of IL-6 were found to be significantly higher with gen- 
otypes containing the G allele in healthy subjects, suggesting a gen- 
etically determined difference in the degree of the IL-6''. In 
multivariate tests, the IL-6 rs 1800795 G/G genotype was positively 
correlated with pain of RA and this may be due to higher IL-6 
responses'". These data have shown that IL-6 rs 1800795 may play 
functional and role in pathogenesis of RA. Although we did not find 
the associations of IL-6 rsl800795 with clinical parameters, further 
study enroll large samples combined with serum determination of 
cytokine network will elucidate the genetic mystery of IL-6 in RA. 

We identified four major haplotype in IL-6 polymorphism: 
GCTGCT, GGCGCT, GGCGGT, GCCGCT. Significant difference 
was observed in haplotype frequencies of GCCGCT between RA 
patients and controls (P = 0.0001, OR = 4.066, 95%CI = 1.891 ~ 
8.746). Haplotype GGCGCT frequencies was found lower in RA than 
controls (P = 0.006, OR = 0.669, 95%CI = 0.501 ~ 0.894). The result 
of haplotype suggests rs 1800796 has potential role in the pathogen- 
esis of RA when in LD with other SNPs. 

A limitation of our study is the weak level of evidence for meta- 
analysis according to Venice criteria. Even though lack of heterogen- 
eity supported by the p-value (0.8028) from egger's test suggests 
absence of critical biases, the accumulated sample size remains small 



for the Asian populations. Although our study was quite large com- 
pared to previously conducted studies, the number of minor alleles 
present in cases in this meta-analysis is rather small as expected from 
the observed MAP. 

We did not find any association between rs 1800796 in the pro- 
moter of IL-6 and RA, which was consistent with a study in a Turkish 
population''. However, this meta-analysis was rather heterogeneous 
with primary studies showing significant results in both 
directions'"'"'^. 

In summary, the data from the present study indicate that a rare 
variant in IL-6 may play a significant role as a genetic susceptibility 
factor for RA in the Chinese population. Further studies with larger 
sample sizes are required to delineate the differences in genetic sus- 
ceptibility to RA in various ancestral groups. Sequencing and func- 
tional characterization of IL-6 may reveal rare variants which may 
influence RA susceptibility in Chinese populations. 

Methods 

Patients and controls. Two hundred and fifty-six Han Chinese patients from 
Shandong province with RA diagnosed according to the 2008 Classification Criteria 
of the American College of I^eumatology and 331 healthy controls from the same 
geographic and ethnic background were included in this study. RA patients were 
recruited from the department of Rheumatology; Provincial Hospital affiliated to 
Shandong University, China. Controls had no history of infectious or chronic 
inflammatory autoimmune diseases and were unrelated to the patients. The study was 
approved by the Local Ethical Committee of the hospital and university, and 
informed consent was obtained from all subjects. The methods were carried out in 
accordance with the approved guidelines. Laboratory parameters were recorded 
including Rheumatoid factor (RF), anti-cyclic citruUinated peptide antibody (aCCP), 
anti-Glucose phosphate isomerase (aGPI), erythrocyte sedimentation rate (ESR) and 
C-reactive protein (CRP). RF was measured usmg N Latex RF Kit (SIEMENS, 
Germany). The level of aCCP and aGPI were determined by enzyme linked 
immunosorbent assay (ELISA) (AESKULISA, Germany). ESR was determined by the 
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Figure 1 | Forest plot of genetic association studies of rsl800795-C and RA in Asian. 



I- 



+ 



0.01 0.1 



1 10 100 
RA Control 



Westergren method via a Monitor- 100 macliine (BALZELLA, Italy) and CRP by 
immune turbidimetric analysis (SIEMENS, Germany). 

SNPs selection. The three promoter region polymorphisms rsl800795, rsl800796 
and rsl800797 in IL-6 had been previously reported as associated with RA in Asian 
population with contradictory results across studies. We expanded the investigation 
in this genetic area by also searching for additional SNPs identified in IL-6 in the 
database of dbSNP (http://www.ncbi.nlm.nih. gov/SNP/). We selected 3 additional 
SNPs with minor allele frequency (MAP) > 1% in Chinese Han population. Along 
with the previously 3 reported SNPs, in total, we included six SNPs in our study; 
rsl800795 (G/C), rsl800796 (G/C), rsl800797 (G/A), rsl524107 (T/C), rs2069840 
(C/G) and rsl474347 (T/G). 

Genotyping. Genomic DNA was isolated by RelaxGene Blood DNA 
System(TIANGEN, China) and the target polymorphisms of IL-6 were amplified by 
polymerase chain reaction(PCR) using TaKaRa Taq Hot Start Version (TaKaRa, 
Japan) and analyzed by HRM methods using LightCycler480( ROCHE, Germany). 
We confirmed our genotyping results by direct sequencing using an ABI3730 
sequencer (ABI, America). 

Statistical analysis. The distribution of genotypes and allele frequencies of RA 
patients for the 6 IL-6 SNPs was compared with control groups using the Pearson's 
chi-square test or Fisher's exact test as appropriate. All p values were two sided and 
using Bonferroni correction for the 6 comparisons significance can be deemed at the p 
= 0.008 level. Hardy- Weinberg equilibrium was examined by chi-square goodness- 
of-fit test using gene frequencies of the healthy individuals. We assessed the LD of the 
SNPs in our sample by using SHEsis (http://analysis.bio-x.cn). We assessed the 
association between the IL-6 variants and RA using a per-allele model and running a 
logistic regression. The distribution of age, disease duration, ESR, CRP, DAS28 and 
autoantibodies with genotypes of each polymorphism with RA was tested using the 
one-way ANOVA test or Kruskal-Wallis test. The association of gender and 
autoantibody status was measured by the Pearson's chi-square test. P values < 0.05 
were deemed significant. LD was investigated using Haploview. 

Meta-analysis. To compute the association between the IL-6 SNPs studied with RA, 
we scrutinized PubMed, Embase, HuGNet and China National Knowledge 
Infrastructure (CNKI) using the keywords 'rheumatoid arthritis', 'RA', and 'IL-6'. RA 
was defined according to the 1987 or 2010 ACR (American College of Rheumatology) 
criteria. The data were last accessed on Dec 30, 2012 and we did not set any language 
restrictions. The eligibility criteria included case- control or cohort studies that 
quantitatively assessed the relationship of an IL-6 gene polymorphism and risk of RA. 
Cases with RA were eligible regardless of age, gender and ethnicity. Studies without 
data on allelic counts in cases and controls and studies focusing on survival or other 
clinical outcome in RA were excluded from the study. In the case of overlapping cases/ 
controls in different studies, we retained only the studies with the largest sample size. 

For each study, we recorded the first author, year of publication, the ethnic group of 
the study population, the method of genotyping, the number of genotyped cases and 
controls, the source of control population, the standard errors (SE) or the 95% 
confidence interval (CI), and the odds ratio (OR). We extracted the calculated ORs in 
each study so as to reflect the same allelic contrast. We performed meta-analysis when 
4 or more studies were available for each SNP. All eligible articles were independently 
screened by two researchers and consensus reached. 

The heterogeneity of the studies was assessed using the Cochran's Q test (con- 
sidered significant for P < 0.10) and was quantified by the F statistic. P ranges 
between 0 to 100% and were considered low, moderate, large and very large for values 
0%-25%, 25%-50%, 50%-75% and >75% respectively. Both fixed effects (Mantel- 
Haenszel) and random effects (DerSimonian and Laird) models were used to com- 
bine the data. Random effects are more conservative incorporating an estimate of the 
between-study variance and thus providing relatively wider 95% confidence intervals 
when heterogeneity exists. Egger's test was conducted with R( www.r-project.org) to 
evaluate the publication bias. The credibility of the associations was evaluated using 
the Venice criteria. Briefly this includes an assessment of the amount of evidence, 
replication and protection from bias that is simplified into three levels of evidence: 



weak, moderate and strong. The analyses were conducted with RevmanS.l package. 
AH p values were two sided. 
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